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D , Abstract 

Binding energy of the 1 state (ortho-positronium) in QED is cal- 
culated using the one-photon exchange Bethe-Salpeter equation in the 
^ I Feynman and Coulomb gauges for different coupling constants a. Cal- 

culations show there is a remarkable difference in values of the binding 
■ energy for different coupling constants in these two gauges. 

^ '. 1 Introduction 

o 
> 



The purpose of this paper is to calculate binding energy of the two-fermion 
1~ system (ortho-positronium) in quantum electrodynamics (QED) using 
the Bethe-Salpeter (BS) equation for the Feynman and Coulomb gauges for 
5^ ! different values of the coupling constants a by a variational method suggested 

recently [1] and see the difference in values of the binding energy for different 
coupling constants in these two gauges. 

The motivation for this study is to look for an acceptable method of in- 
vestigation of bound states in systems like quantum chromodynamics (QCD) 
in which the coupling constant ~ 0.1 -i-0.5 is not too small as in the QED 
a = ^ = 0.0073. In general, it is generally agreed that the BS equation in 
the one-photon or one-gluon approximation appears as a quite acceptable in- 
strument for study of a bound state problem if the coupling constant is small 
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enough. In principle, it is hoped that this approximation gives the main con- 
tribution to binding energy. However, the BS equation in gauge theories is 
not gauge invariant in this approximation and in QCD the couphng constant 
is not so small as in QED. Therefore, we should know what difference in 
binding energies arises for different gauges and different coupling constants. 

Most computations of electromagnetic bound states use the Coulomb 
gauge and the Breit potential with relativistic corrections (see, for example, 
[2, 3, 4, 5]). This approach is most compatible with experimental data. 

Long time ago an attempt was made [6] to consider gauge properties of 
the BS-cquation for the two-fermion electromagnetic bound state for different 
covariant and axial gauges. Naturally, it was found that different gauges give 
different results in the highest orders in electromagnetic coupling constant 
a. 



2 Lagrangian. 



Wc perform all calculations in the Euclidean metrics. The QED lagrangian 
describing electrons and photons looks like 



L = -^F^^(x) + {iIj{x){p + eA{x) - m)ij{x)), 
F„^{x) = d^A^{x) - d^A^{x). 

The electron propagator has the standard form 



S{x-x') = (t 



dp e'^^^ ^'^ 
{2%)'^ m — ip 

The photon propagator is defined by the formula 

dk 



D^,{x) = (T[A^(x)A,(0)]) 



(27r) 



jD^^{k)e 



ikx 



where in the Feynman gauge 
5, 



D,.{k)^^-^, D,^{y)^S,.D{y), D{y) - 



(1) 



(2) 



(3) 



(4) 
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and in the Coulomb gauge 



'W 



J_ 

fe2 



< 



(5) 



j_ 

k2 



2.1 States of electron-positron system 

Let S be spin and L be orbital momenta. Then the total parity of electron- 
positron system is defined as P = (—1)^^''". For para-positronium S — 
0, = 0^ and for ortho-positronium .S* = 1, J^' = 1~ . 
Quantum numbers of electron-positron currents 



can be determined in the nonrelativistic representations of electron and positron 
wave functions 



Quantum numbers of all possible electron-positron currents are represented 



(*(-k)Oj*(k)) - {v{-]<.)Oju{k))b{-k)a{k) 



(6) 



^(k) 



u(k)a(k) +v(k)b+{k), 
Il(k)a+(k) +U(k)6(k), 



where 




in Table 1. 
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Table 1. Quantum numbers of relativistic currents 
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1+ 
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-1 


0- 
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+1 


1+ 
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(*75*) = 


^ 1 











-1 


0- 



One can see that only vector and tensor currents have the ortho-positronium 
quantum numbers. In other words the ortho-positronium should be described 
by a mixture of vector V and tensor T relativistic currents. 

2.2 One photon exchange and ortho-positronium cur- 
rents. 

In order to extract the currents having ortho-positronium quantum numbers 
in the one photon exchange approximation it is convenient to use the method 
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of generating functional. We have 

(7) 



where 

^2 



Let us introduce the new variables 

.y y 

X2^X-- 

The Firz transformations should be used to extract the vector and tensor 
currents 



pa 



Yl {Oji)aa{Oj^)pi3^TT[Oj,i'y^OjJj^]D^^{y) 

Jl,J2 



3 Feynman gauge 

In the Feynman gauge (4) we have 



Thus, in the Feynmam gauge the ortho-positronium is described by the vector 
current only. 

For the vector-vector part of the one-photon exchange contribution one 
can get 



W2 =^ ^ jj dxdyJi{x,y)Vij{y)Jj{x,-y) 



e2 



jj dxdyJj{x,y)D{y)Jj{x,-y), 
Mx, y) = (x + I) 7,^1^ - I)) = (^) {x) 

Px= - (^dx - dx^ 

Let us introduce an orthonormal system of functions 

f Jdy U*Q^{y)UQ^{y) = 5q,,q„ 

{UQ{y)} = {Un.im{y)} = \ ^ ^ ^^^^^ ^ ^ (10) 

E UQ{y,)U^{y2)^5{y,-y2), 
K Q 

As long as for the ortho-positronium state Q = we restrict ourselves to the 
function Uo{y) = U{y) with normalization {UU) = 1. 
Let us perform the following transformations: 

W2 = ^ J dx JJ dyidy2Ji{x,yi)^D{yi)5{yi-y2)^D{y2)Jj{x,-y2) 

dx Jj{x)Jj{x), 



e 

T 



J^{x) = |rfyy^C/(y)J,(x,y) = (*(a;)y(p,)7,*(x)), (11) 
where the vertex is 

V{Px) = Jdy .jD^)U{y)e^^^% (12) 

The generating functional Z containing the ortho-positronium vector current 
can be transformed as 



with 



2 

-^-TYlV-iiBiSV^jBjS] dxidx2Bi{xi)Uij{xi - X2)Bj{x2) 



The polarization operator is defined by 



Uij{xi - X2) = - y'^^ 



V (P^,) -fiS{Xi - X2)V (p^,) -fjS{x2 - Xi) 



dq 



where 



dk 

2 J (2^ 



V'^ik) Tr 



liS{k + -]jjS[k- 



dk 

2 y (2^ 



V{k)- 



Tr 



7j + ik + z|j 7j (^m + ik 



''2 



m' 



dk 

2 i (2^ 
2e^ ■ 



A:2 - |k2 



W (m2 + (A;+|)') (^m2+(/c-|)' 



•5,,- (14) 



We choose the frame where q — (0, iM) . 

According to [1], we divide the polarization operator (14) in two parts 



flijiq)^ Uo{q)+ni{q) 



where 



dk 



J f27r 4 ^ ' 



V\k) 



m- 



+ {k+l) 
(ik^ + kl) 



m- 



(k'^ + kl + - + M^kj 
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is responsible for a continuous spectrum, and 



2e' 



I 



dkdk. 



(k2 + A;| + m2-Mi)' + M2A;|' 



is responsible for a bound state. 
The vertex (12) looks like 



V{k)^V{k,h)^ dy dy^ 



*(y,?/4)e 



jyk+ifc4{/4 



(27r)2(y2 + y|) 



with U{y) — ^jD{y) ^(y,y4), where ^(y, 2/4) is a wave function. 

Next we should extract in the kernel of (13) the part that is responsible 
for the ortho-positronium bound state 



{B[I - Il]B) = {B[I - Ho - n,]i3) ^ B 



I - 



J-Hn 



B 



Finally the binding energy is defined by the equation: 



1 — max 



max 



{VIliV) 



V 



max 



(^ [D - dUqD] ^) ■ 

(15) 



3.1 Nonrelativistic limit 

Let us consider the nonrelativistic limit in the equation (15) which takes 
place for small coupling constant. We should take 



*(y,y4) = *(y) 



with 
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Equation (15) becomes 



For the polarization operator one can get 

Uj(q) = 2e^ / - — -fy2 k,A;4 ^ — 

(2^)' {k^ + kl + m^-^Y + M^kl 



. f dkdk^ \/-(kJ ^ ^2^ /• 

i (2ttY (\s? + me? + Am^k'i J 



(27r)4 (k2 + mey + Am^kj J {27rf + me 

In order to get the nonrelativic Schrodinger equation, the term A;| in the 
denominator should be neglected (for details see [7]). 
Finally, we have the equation 

jrW^(y)rAT^*(y') 

1 = am max ■ 



which is nothing else but the nonrelativistic Schrodinger equation. 
In the non-relativistic case 

*(y) = e-"^ (17) 

and we get the well known result 

1 = 2am max . . „ =^ e = — m. 

» (a + y/mey 4 

3.2 Variation calculations 

Let us come back to equation (15). For the lowest state of orthopositronium 
the test function is chosen in the form 

*(y,|/4)=e-«V^^, (18) 

where a and 6 are variational parameters. This function is the closest to the 
nonrelativistic wave function (see [7]). 
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Then the vertex function is 

V{k) = JdyD{yMy,y,)e-^'y = Jj^^^e-^V^^.-^^^^^^^ 

a6 + v^^^^ITWT^ 



and 



[k2 + s2 + a2(l + b)+ 2aVa% + k% + s2] Va^^ + k^fo + s2 
R{k,s;a,b), k = \k\, s = |/c4|; (19) 



J (27r) 



(27r)2t/2 4^2(1 + ^) 

The positronium binding energy e = 2m — M is defined by the equation 



^ 4m2y a,b TT^y y + + A)2 + 4(1 - A)s2- 

A = l--— = -(1--^ , M = 2m-e. (20) 
4m^ m \ Amy 

Since of A is small, it is convenient to introduce the new variables 
k VA k, As, a ^ VA a, \/A b. 

We get 

H{k,s,A;a,b) 

VAab + + \s?b + s2 



[k2 + A s2 + a2(l + Ab) + 2a^A{a?b + \i% + s^)] ^^2^+k26T^ 
and our equation takes the form 

^y-^ Am^J aTiT^l I (P + As2 + i)2 + 4(i_A)s2 • 

(21) 



Preliminary calculations have shown that the parameter a is very close to 
one and the parameter b is very small for all coupling constants a < 0.5, so 
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that we can put a = 1 and 6 = in the hmits of our calculation accuracy. 
Thus, the test function practically coincides with the non-relativistic wave 
function *(r) = e"^ ^ 

In the case a—1, b — the function H is 

H(k, s, A; 1, 0) = 

Finally, the equation defining binding energy of the ortho-positronium 
looks like 

'^7KV-2)^^J/^J/'- (fc2 + A.2 + l)2 + 4(l-A>2 -^''^ 

One should stress that the difference between relativistic and non-relativistic 
cases is defined by the term — A s"^ in the denominator of the fermion 
loop. 

The results of numerical calculations are shown in Table 2. 

For semiquantitative calculations one can use the approximation 



^'J J (F + A5^ + l)^ + 4(l-A).= 2^1 + 9x/A + 6A 

SO that the equation become the form 



1 _ 2 - A ^ _ 4^A(l + 9yA + 6A) 

^ 4^A(1 + 9VA + 6A) 



or a 



2- A 



This formula gives semi-quantitative dependence of the binding energy A on 
the coupling constant a. 

4 Coulomb gauge 

Now let us consider the Coulomb gauge. In this case the Firz transformations 
lead to 
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= Yl iC>Ji)aa{0j^)pi3^T:^[0j,i^f,0jJ-f^]D^^{y) 

Jl,J2 



-{7j)aa'DYr{y){7j')pP - {il0lj)aa'DL{y){ijQjj')pi3. 
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where 



1 

4 



k2 



1 



4 k2 







1 ^ 1% 






1 4k2' 






1 . 


k2 




1 4k2' 



(23) 



It means that in the Coulomb gauge the ortho-positronium is described by 
a mixture of vector V and tensor T relativistic currents. 

In what follows we neglect the term There are two reasons to do 

it. First, usually in the generally accepted approaches these terms are not 
considered at all. Second, we did not have courage to perform these cumber- 
some calculations although they can be done in case of emergency. Thus we 
have 

1 5iA 



and 



4k2 



167r 



The vector and tensor currents are 



The one-photon exchange term containing vector and tensor currents looks 
like 



Px^ - { dx - d, 



^2 = ^11 dxdy 



J^i^, y) ■ '^luiy) + Jjii^, y) ■ '^luiy) ■ Ju i^, -y) 



- • 5(yi - y2) ■ 



yii 



y2| 



+ 



S{yi - y2) 



yi| 



y2| 
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In the Coulomb gauge we should introduce an orthonormal system in the 
X G R'^ space only: 



{UQ{y)} = {Unimiy)} 



i I dy U*Q^{y)UQ,{y) = 5q„q,, 

E UQ{y,)U*Q{y2) = 5{yi-y2), 
k Q 



For ortho-positronium Q — and Uo{y) — U{y) with normalization 
{UU) = 1. 
We have 



Here 



J^{x) = (^{x)v(p,)rT^{x)) , Tt = i7o7 



with 



dy ^e-^^ 



The generating functional takes the form 



where 



DB = DBvDBt, (BB) = (ByBy) + (BtBt) 
B = VTvBv + VTtBt. 

The integration over electron variables gives 

J I D*D^e(*^"*)+^(*^*) = exp |Tr In 1 + ^BS 

K, e"t // dxidx2TrBS{xi-X2)BS{x2-xi)] = JI (^{B(p)n(p2)B(p)) 



(24) 
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and the generating functional becomes the form 

Z = ||DBe-5(BB)+^(BnB) (25) 

Here 

i^p ) - \^ n('^^)(p2), n(^^)(j92) J ' ^^^^ - \ Bt(p) J ' 
n(^i^2)(p2) ^ ^ / ^ ^ ^ ^ '^^ 5^ (26) 



7;5;^^^^)(fc,p) = -^Tr 



Direct calculations give 



rji,i f m + A: + I j Tj^j ( m + A; - | 



(27) 



^0 



2,2 , 
Av4 i- 3 







"'4 "o" 



imM 

M2 



-imM + ^ 





The matrix Tq is responsible for a continuous spectrum and should be re- 
moved. The matrix Tj is responsible for bound states and can be represented 
as 

2 



m + 



M 



The eigenvalues are 



1 

—i 



A 



(1 z) / M 
^ ^ + im- — 



± 



(28) 



For the eigenvalue A = (m — —171^(1 — = (^2^j <S the 
bound state can exist but this case requires a very large coupling constant a 
to provide the condition 

TT TT \2mJ 
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This condition cannot be realized (see [1]). 

The most favorable situation for the existence of the positronium takes 
place for 



The quadratic form in the representation (25) can be written in the form 
(BB) - (BHjB) = {W+[l - Il+]W+) + {Wl[l - n+]V^_), 

W±=-^(By±iBT), 



n+ = a 



dk V\k){m±fy 



The combination = ;^(By+iBT) with 11+ describes the desirable ortho- 
positronium 1~ state. The Dirac spinor of this state equals 

«+=(l + 7o)«=(^Q oj(^_^)«^ = 2(^jjw 

according to a standard approach (see, for example, [2]). 
Finally, the partition function reads 

The ortho-positronium bound state equation looks like 

1 = n+ (29) 

4.1 Numerical calculations 

For numerical calculations let us choose the test function in the form 

u{y) = ^ = \\u\\' = (uu) = ^ 
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so that 

V'(k) = J- /<iy ^e'(M =^[iy L^e'^^' = J^^^, 
11(711 J / y V |y| + m^'a^ 

In our calculations we have only one variational parameter a. 

After some calculations equation (29) can be represented in the form 

l = maxn^ = max^(^l + — j —j ^ ^ 

where 



2V2 



[■\/2c2 + 62 + 5V4c2 + 62 + ^2c2 + 62 _ 6V4c2 + 62] 

The results of the numerical calculations are given in Table 2. The accuracy 
of these calculations is about 1 -j- 2 % (see [7]). 
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Table 2. Binding energy e (eV) of the state 1 as a function of the couphng 

constant a. for different gauges. 



a 


0.0005 


0.001 


^ = 0.0073 


0.01 


0.1 


0.3 


0.5 


Feynman 


0.032 


0.126 


6.47 


12.0 


893 


5 700 


12 600 


Coulomb 


0.032 


0.127 


6.8 


12.8 


1 270 


10 800 


27 800 


Schrodinger 


0.032 


0.127 


6.8 


12.8 


1 280 


11 500 


31 900 



5 Conclusion 

In conclusion, let us formulate our results. 

• The Feynman and Coulomb gauges give coinciding results for very small 
coupling constants a < O.laQED- 

• For aQED = yI? difference of binding energies for the Feynman and 
Coulomb gauges is of an order of ~ 5%; 

• For a ~ 0.1 -i- 0.5 difference in the binding energies is of an order of 
~ 100%. 

• Calculations in the Bethe-Salpeter equation with the Coulomb gauge 
and in the non-relativistic Schrodinger equation coincides up to o; < 
0.1. 

Thus, one can conclude that in gauge theories hke QCD, where the couphng 

constant is not too small, the Bethc-Salpeter equation in the one-"gluon" 
exchange approximation gives quite different numbers for different gauges 
and, therefore, it is not a good mathematical instrument for calculation of 
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binding energies of bound states. One can say that the gauge invariance is 
broken in the Bethe-Salpeter equation with any fixed kerneL An alternative 
way is to recognize that there exists a preferred gauge, namely the Coulomb 
gauge. This idea is not new (see [8] and other references there). 

Besides, in real QCD we have an additional difficulty: the formation of 
mesonic bound states takes place at large distances where confinement plays 
the main role and we should know the explicit form of quark and gluon 
propagators in the confinement region. 

The author is grateful to E.A.Kuraev and V.N.Pervushin for useful and 
stimulating discussions. 
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